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Abstract 



Fluids of grade n are continuous media in dynamic changes of phases avoiding 
the surfaces of discontinuity and representing the capillary layers in liquid-vapour 
interfaces. We recall the thermodynamic form of the equation of motion for inviscid 
fluids of grade n [1]. First integrals and theorems of circulation are deduced. A 
c/2 • general classification of flows is proposed. 
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O ■ 1 Introduction 



In continuum mechanics, media of first gradient order cannot represent fluids 
with strong density variations. Material surfaces need their own characteristic 
behaviour and energy properties [2]. 

Fluids of grade n (n is any integer greater than 1 ) are continuous media with 
an internal energy per unit mass e which is a function of the entropy s, the 
density p and spatial gradients up to the n — 1 order: 

e = e (s, grad s, (grad)*™ -1 ^, p, gradp, (grad)^ _1 ^p 
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D.J. Korteweg pointed out the advantage of fluids of grade n with respect to 
the density and removed the discontinuity surfaces representing the capillary 
layers in liquid- vapour interfaces [3]. 

We improve the model accuracy by taking into account the successive gra- 
dients of density. Recently, this approach was used in the case of dynamic 
changes of phase [4]. Till now, the thermodynamic part of the fluid has been 
neglected. Due to the fact it is not possible to consider a virtual displacement 
of temperature, it is necessary to use the entropy by the way of the internal 
energy. 

The equation of motion with both strong gradients of density and entropy 
is obtained. One can accord importance to the fact that the thermodynamic 
form of the equation of motion is independent of the order n. Flows with 
strong variations of temperature such as those associated with combustion 
phenomena or non-isothermal interfaces may be considered. 

The conservative flows of classical perfect fluids are only a mathematical 
model. It is also the case for fluids of grade n. The study of their structure 
is necessary. The conclusions for fluids of grade n are of the same kind as for 
compressible fluid [5]: it is possible to obtain a Clebsh transformation of the 
motion equations {potential equations) and many invariance properties and 
first integrals. The potential equations may be used to classify motions in the 
same manner as in the case of perfect fluids [6,7]. 

The simplest case is for fluids of grad 2 (or thermocapillary fluids) : the inter- 
nal energy is a function only of the density, the entropy and their first spatial 
derivatives 

e = e(s, grad s, p, grad p) 

and n > 2 adds little except for complications in the momentum balance 
equation. If we consider flows through liquid-vapour interfaces and use conve- 
nient physical units, the equations of thermocapillary fluids are able to study 
non- isothermal motions and yield a possible interpretation of film boiling phe- 
nomena [8,9]. 



2 Thermodynamic form of the equation of motion for inviscid flu- 
ids of grade n 

The principle of virtual works allows to obtain the equation of motion. In the 
case of a perfect fluid motion, it is written by means of the Hamilton principle 

[5,6,10]. ' ' ' 

The variations of particle motions are deduced from families of virtual motions 
of the fluid in the form [5] 

X = *(x,£, a) 
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where a is a parameter defined in a neighbourhood of zero. The real motion 
corresponds to a = 0. 

A particle is represented in Lagrange coordinates by X (X 1 , X 2 , X 3 ) as a po- 
sition in a reference space V Q . At time t its position is given in V t by the 
Eulerian representation x (a; 1 , x 2 , x 3 ) (see Appendix). The virtual displace- 
ments associated with a variation of the real motion can be written 

5X = —(x,t,a) 
da 

Q = 

The variation is dual with Serrin's [5]: the two variations are mathematically 
equivalent [11]. The Lagrangian of a fluid of gradn is defined as 



,(Iw- e -n) 



where V denotes the velocity, Q denotes the extraneous force potential defined 
in T>t and * is the transposition in T> t such that V*V is the scalar product of 
V by V. Between times t\ and time t 2 , the Hamiltonian action is 



rt 2 

L dvdt 

Iti JD t 



Jt! Jl 

The density verifies the balance equation 

pdetF = Po (X), 

where p a is defined in V Q and F is the deformation gradient. We deduce 

^ + p div V = (1) 

The entropy variation is the sum of a variation associated with the virtual 
motion and another one, 5±s, related to the particle [9] 

where Sis and 5X are independent. To obtain the equation of motion, we con- 
sider the conservative case corresponding to S±s = 0. 

Classical methods of variational calculus provide the Hamiltonian action vari- 
ation 

5a = 

A calculus as in [1] yields 

A 1 

5X dv Q dt, 



5a = f f p 
Ju Jv t 



d 

— iy*F) — #grad* s — grad* m 

CLL 
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where 



m = -V*V -h-Sl 
2 

h = e + P - 
P 

n-l 



P = p 2 e' p + pj:(-l) k div^\, 

e 



de 



8 (grad (fe) p) 

4 + 1 E(-i) fe div«L de 
p k=i v ( 



(2) 



1 dp Q 



5~K and 



(grades) 

We have taken into account the relation bp = pdiv ((5X) + 

det F oX 

5V = — F — -4 — -; div^ and grad^ are the divergence and the gradient op- 
dt 

erators in V t reiterated A;-times; div Q and grad D are the divergence and the 
gradient operators in V . 

Obviously, p, 9 and h have respectively the dimension of a pressure, a tempera- 
ture and a specific enthalpy. We call these quantities the pressure, temperature 
and enthalpy of the fluid of grade n . 

The Hamilton principle yields: 

For any displacement 5X null on the edge of V Q , 5a = . 

We get the equation of motion 

d(V*F) 



dt 



— 9 grad* s — grad* m = 



dV d(V* F) 
Let us note that (T* + V*— )F = ^— — -, we obtain 

ax dt 

r = 9 grad s - grad(/i + SI) 



(3) 



In case of inviscid thermocapillary fluids the dissipative function is null and 
the equations of motion and energy imply [12]: 



. ds 

p9 — + div q — r = 

dt 



(4) 



where q is the heat flux vector and r the heat supply. This result is extended 
to the fluids of grade n [13]. 

A constitutive equation must be added to these equations, which yields the 
behaviour of the heat flux vector q. For example, the Fourier law is 



q=— /cgrad# and k>0 



(5) 
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The heat supply is assumed to be given in similar way than the extraneous 
force potential. From the heat conduction inequality, we obtain the Clausius- 
Duhem inequality and the fluids of grade n are compatible with the second 
law of thermodynamics [14,15] 



ds /q 



r 



Consequently, the equations of motions of inviscid fluids of grade n verify 



r = 9 grad s — grad(/i + Q) 
dp 



dt 



+ pdivV = 



n ^3 

p 9 — + div q — r 

dt 







(6) 



q=— fcgrad(9 and k>0 



The thermodynamic form of the motion equation of inviscid fluids is equivalent 
to the classic balance equation. The thermodynamic form is the same for all 
fluids (independently of the order n). Nevertheless, the complexity of balance 
equations increases with n. In fact, only the expressions for the temperature 
and the enthalpy become more complex and we have to modify the definitions 
of these quantities^. 



3 Thermodynamic hypothesis on the motion 



If we study the invariance properties of the motions, we must assume that the 
heat supply allows a behaviour expressing that a thermodynamic quantity T 
has a zero material derivative or - with a stronger assumption - is constant 
everywhere in the flow. Consequently equations (@J and ([5]) are replaced by 
one of the two conditions 

jrri 

— = or T = C ste (7) 

dt v ; 

For example, we will consider the case of an isentropic motion (T — s). Then, 

dt 

For the case of an isothermal motion (T = 8). Then, 

9 = C ste 



1 The definitions are given by system (|2|). 
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They are two limit cases: one is associated with fast motions, the other with 
slow motionj 2 "! Then, the fluid motions satisfy the system consisting of equa- 
tions (|U), © and 0. 

We consider the more general case when T is a differentiable function of 9 
dT 

and s such that -— ^ 0. Consequently, we denote by H(T, s) a differentiable 

OH 

function of the variables T and s such that — — = 6(T, s). 

OS 

The function H is defined, to an unknown function of T, by the differential 
form 

dH = 6 ds — u dT 

The two variables u and T can be used instead of 9 and s to study the general 
case by transforming Eq. (J3j) following the method of Lemma 1. Let us denote 
g = h — H , Eq. (ED can be written: 

T = u grad T - grad (g + Q) (8) 

Let us note that if T = 9, g is the free enthalpy (g = h — 9s) and u = —s, and 
the motions satisfy System ( (DO),©, (jHJ)). 

First integrals are obtained by the research of invariant scalars as the state 
variables of the fluid (the density, the entropy, the velocity, with a zero material 
derivative) and expression ([8]) allows us to make a single classification for the 
fluids. 

The momentum equation expresses the balance of forces (inertial forces, body 
forces and stresses). 

Let us note that System ((JI]),©,©) is not in balance form, but expresses the 

dT 

fluid evolution: the relation — — = implies T is constant along trajectory and 

dt 

the equation of continuity — = — pdivV yields the density variation. It will 

dt 

be shown that equations (j3J) or (jBJ) describe the evolution of the velocity field. 



2 The Euler equation represents the equation of motion for classical inviscid flu- 
ids. Such motions can be obtained with non-negligible heat flux and heat supply. 
They correspond to physical situations. As a matter of fact, in the non-dimensional 
equation of motion for viscous fluids [16,17], it is often justified to neglect the vis- 
cous terms when, with respect to other physical quantities, the Reynolds number 
is large. Nevertheless, for large temperature gradients and with a Prandlt number 
of order one, it is not possible to remove the terms of thermal conductivity. Such a 
case arises when we study the "thermal" boundary layer [18,19] and slow motions 
of natural convection in fluid with gravity and heat sources associated with large 
differences of temperature [20]. 
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4 Material and convective derivatives 



The material derivative of any tensorial quantity expressed in the physical 
space is the derivative with respect to time when one follows the particle 
along the motion. 

The tensorial quantity has an image in the reference space. The image in the 
physical space of the derivative with respect to time of the tensorial quantity 
in the reference space is the value of the convective derivative. The convective 
derivative is denoted by d c (see Appendix). 

The material and the convective derivatives of a scalar are equal. A tensor is 
convected by the flow when its convective derivative is null. 

If we transpose Eq. (151) , we get the convective derivative of the velocity covector 
V* in the form 

d c Y* = grad* (- V 2 - g - Vt \ + u grad* T. (9) 

Eq. ([9]) is equivalent to Eq. ([S]) and yields the Helmholtz- Kelvin theorem in its 
most general form: 

Let us consider an arbitrary fluid curve C such that I{t) = J V*dx is the 

circulation of the velocity field. The circulation of the convective derivative of 
the velocity field is 

*L = [ d c V*dx = { d(-V 2 -g-tt) + { udT 
dt Jc Jc \2 J Jc 

If C is a closed curve, the integral J d ( — V 2 — g — Q^j is null. Moreover, if 

C is located on a T-surface (surfaces T = C ste are fluid surfaces and this 
assumption is compatible with the fact that C is a fluid curve), the second 
integral is null and I is constant along the motion. 

When the entropy is constant all over the fluid (T = s = C ste ), the so-called 
Kelvin theorem corresponds to the special case of homentropic motions: the 
circulation of the velocity field on any closed curve is constant. 

Due to Noether's theorem, it is known that any conservation law can be rep- 
resented by a group of invariance. It has been shown the conservation law ex- 
pressed by Helmholtz-Kelvin theorems associated with isentropic fluid curves 
corresponds to the group of the permutations of the particles of equal entropy 
[21, 22]. This group keeps the equations of motion invariant for both a classical 
perfect fluid and a perfect fluid of grade n. We propose to define a general per- 
fect fluid as a fluid invariant by this group or as a continuous medium whose 
motions are submitted to Helmholtz-Kelvin theorems. 
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5 Potential equations 



Let p and ifj be two scalar fields such that the material derivatives verify 

—— — - V — q — it and — — = u. 
dt 2 y eft 

Equation (jUj) is equivalent to 

d c V* = 4(grad* 9? + ^ grad *T) . 

The covector V* and the covector W* = grad* tp + ip grad* T have the same 
convective derivative. The difference (V — W)* is a covector convected by the 
flow, which leads to the system 

V = grad p + tp grad T + C, 

(10) 

dT dip 1 2 dip 

^ = ' -dJ=2 V ~ g -^ Tt= u > 4C =0 - 

System (([I]),©,®) is equivalent to system ((fTl),(JT0l)). 

In this case, the generalized Helmholtz-Kelvin theorem is reduced to the cir- 
culation of C along any fluid curve is constant. 

Due to the Stokes formulae, the vector U = - rot C is convected by the flow 

P 

(d c U = 0) i.e. : 

dt dx ^ 
Due to the fact that C = 7 grad is, where 7 and v are two scalars convected 
by the flow and grad p, grad^, grad v are three independent vectors (see 
Lemma 2 and Theorem 3), the system (TTU]) leads to the system of potential 
equations 

V = grad p + if) grad T + 7 grad v, 

(12) 

? = 0, ^ = -V 2 -,-fi, d 4 = u, 5=0, 5 = 0. 
dt dt 2 ^ dt dt dt 

In the general case, the velocity field is associated with five potentials; T is 
one of therrJiD- 



3 The velocity field depends also on the density which is not considered as a po- 
tential. 
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6 Classification of motions 



As for T, the potentials 7 and v have a zero material derivative. They are first 
integrals of the motion [23]. 

These potentials have constant values on each trajectory; when they are con- 
stant in a fluid domain, the values remain constant in time-space. 
In the domain convected by the flow, the fluid velocity is expressed by fewer 
potentials characterizing the kinematic properties of the motion. In an adjoin- 
ing fluid domain, the motion may be different. Then, the flow is separated into 
different domains that move, remain adjacent and keep their own character- 
istics. 

Due to the absence of viscosity, in hydrodynamics of incompressible perfect 
fluids the flow is divided into vortical and irrotational domains that do not 
mix. The following study shows that, for all inviscid fluids of any grade (com- 
pressible or not), the vortical domain may be divided into several parts of 
different kinds 4 1 

First integrals are tensorial quantities of different nature. For a non-scalar 
tensor, the first integral is a tensor convected by the flow (or with a zero 
convective derivative) i.e. has a constant tensorial image in a reference space. 



6.1 Oligotropic motions 



When the covector C* is zero (i.e. 7 = or v = C ste . ), the velocity field has 
a particular kinematic type. We have the following properties: 

(a) The velocity field can be expressed by the use of three potentials and 
the equations of motion are 



V = grad ip + ip grad T, 



dt ' dt 



- V — q — si, — — 

2 J ' dt 



u. 



(13) 



(b) From System ([IB"]) we get 

rot V = grad ip A grad T, 

then 

grad* T rot V = 



(14) 



System (|12|) assumes the flow to be without discontinuity. A shock wave trans- 
forms the nature of flow [11, 24]. 
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Consequently, the vortex lines are located on the surfaces T = C ste . 

Property (j!4p is equivalent to System (ITS]) . The motions, we call oligotropic 
motions, are intermediate between irrotational and general motions. In System 
(JED, the expression of the velocity is independent of time and the motion 
remains oligotropic at any moment. 



6.2 Homentropic motions 



The thermodynamic condition is T = C (which is stronger than — = 0) 

dt 

and the velocity field is of another particular kind. When T = s , the motion 
is not only isentropic but also homentropic. We have the following properties: 

(a) The velocity field can be expressed with the use of three potentials and 
the equations of motion are 



V = grad <p + ip grad u, 



G?7 



dv 



dT dip 1 9 „ 

l *=°- i=2 V -*- n - *=°' *=°- 



In fact the Clebsch representation [16] is practically worthless because the 
curves 7 = C ste and v = C ste are the vortical lines whose complexity is well 
known [25]. It is easier to write the motion equations in an equivalent form: 



V = grad (p + C, 

< 

I dt ' dt 2 y 



(b) Vector - rot V is a convected vector 
P 

d c (j rotV ) = ( 15 ) 

When a thermodynamic potential T is constant in the flow, the Cauchy the- 
orem for homentropic motions is generalized and Eq. ([15!) is characteristic of 
these motions. 
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Fig. 1. For perfect fluid motions, the fluid domains are distributed in four parts. 
Convected by the flow, the domains become deformed but do not mix. 

6.3 Irrotational motion 



If an oligotropic motion is a motion where T is constant, the velocity field 
derives from only one potential 



V = grad (p, 



dtp 
~dt 



At every time, rot V = GO- This is the case of irrotational motions: the 
velocity whose the evolution corresponds to the equation of Bernoulli is the 
gradient of a single potential. 



6.4 Permanent motions 



For a permanent motion, Equation (|S]) may be written in a generalized Crocco- 
Vaszonyi form 

rot V AV = m grad T — grad G (16) 
where G = - V 2 + q + Q. 

If T — s, g is the specific enthalpy and G the Lamb function [16] (surfaces 

5 When at one moment rot V = 0, this property is not necessarily verified at every 
time as in hydrodynamics of compressible fluids. 
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where G = C ste are Lamb surfaces). Consequently, 

dT „ dG 



dt 







dt 







G is a scalar first integral of the motion and we use G for function v. Potential 
equations for general permanent flows are 



dt 



V = grad <p 1 + ip 1 grad T + r y l grad G, 



1 , , ^ dib dT d'y-, dG 

-V 2 -g-Q, -i- = u, — = 0, -^ = 0, — =0. 
2 y 1 dt 



dT n dG n 
Let us note: — — = and — - — = 
dt dt 



dt 
dt 



dt dt 
= —G; if ip — ip ± + tG, (p is 



independent of time. Let us denote n f l — 7 + 1 and ip 1 = ip, the equations of 
motion can be written with potentials independent of time 



V = grad <p + ijj grad T + 7 grad G, 



dib dT d'-f dG 

— = u, — — 0, — = -1, — = 0. 
at dt dt dt 



Let us note: C = (7 + 1) grad G and d c C* = 0; the behaviour of <p 
is a consequence of the previous equations. 



.dip 
~dt 



(17) 



V 2 ) 



The two scalars T and G characterize the motion. The classification of motions 
is in accordance with the respective positions of the surfaces T = C ste and 
G = C ste (when T = s they are the isentropic and the Lamb surfaces). We 
have different cases 

(a) General case: T and G are independent (or gradT A gradG 7^ 0). 

The T-surfaces (isentropic surfaces) and the G-surfaces (Lamb surfaces) are 
different. 

(b) T is not constant and G = G(T): the T-surfaces are G-surfaces. 
Replacing if) par ip + 7 G' T we get 

V = grad (p + ip grad T, 



(18) 



dip nl(rr , dT 

—— = u — G (T), — — = 0. 
dt v ; ' dt 



These equations characterize permanent oligotropic motions [6]. 
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General mot ions : G and T independent 



Motions G 




( isoenerget ic motions ) 








Motions T = T 









Fig. 2. Classification of permanent motions: the different domains do not mix along 
the motion of the fluid. 

(c) G is constant and T is not constant. 
This is a particular case of (b): (-^ = u). 

(d) T is constant, G is not constant. The equations are 

V = grad (p + 7 grad G, 



(19) 



dl dG m m 

dt ' dt 



(e) T and G are constant. 

We obtain the irrotational motions 

V = grad^, G = G , T = T Q 

Let us note: All the irrotational motions are not represented by case (e) (such 
is the case of Hamel's motions [26]). In fact, relation (TlBl) points out that if the 
motion is irrotational and if one of the mappings T or G is constant, the other 
mapping is also constant. If neither G nor T are constant, G must depend on 
T (case (b)): an irrotational motion is oligotropic. 

Irrotational motions do not correspond to particular values of the first integrals 
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(or values convected by the flow). In Fig. 2, we represent the flow of a fluid 
with different domains in T> t . The domains are convected by the flow and do 
not overlap. If each domain has its own color, it keeps the color all along the 
motion. 



7 Appendix 

The motion of the medium consists of the t-dependent C 2 -diffeomorphism 

<Pt [5] 

X e T> — > x6D, 



We deduce — — = ^— F and ^— — 
dt ox dt 



—F 



_ x dV_ 
<9x 



Let us write T*(T> t ) for the cotangent fiber bundle of T>t and T*(T> t ) for the 
cotangent linear space to T> t at x; then, 



L(x,t)ET*(D t 



represents a field of differential forms on T> t . 

Let us write T*(T> ) the cotangent fiber bundle of T> and T^_(T> ) the cotangent 
linear space to T> at X. The mapping (p\ is induced by ip t for the form fields. 
The convective derivation d c of a form field L is deduced from the diagram 



L E T*(V t 



i* 



LF E T*{V ) 



d r 



dL L &V 

dt dx 



d_ 

dt 



(20) 



Let us note that — — h L ^— is the Lie derivative of L with respect to the 
, at ax 



velocity field V 

The velocity field V is the infinitesimal transformation of the one-parameter 
group of transformation ip t [27]. 



Generally, for any tensor, the convective derivative is the Lie derivative with 
respect to the velocity field. 
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Consequences : T> t is assumed to be an Euclidian space; let b be a scalar field 
on T> t . We define two form-fields as V* and (grad &)*. From (I20I) we deduce 

*< v *>= r,+ sG v2 

where T is the acceleration vector. Taking into account that 

4(grad£>)* = (grad^)* and d^LF" 1 ) — ^ F" 1 , 
Eq. © yields 

4V* = grad* Qv 2 - £/ - ^ + 6>grad*s. (21) 

Lemma 1 Let 9 and s be two differential scalars fields in T> t . 
There exist three scalar fields a, (3, 7 such that 

9 grad*s = a grad*/3 + grad*7 (i) 

is equivalent to: there exists a function ip such that 

P = <p{9,8) (ii) 



Proof : 

- If grad 6 A grad s = then 9 = $'(s), 9 grad s = grad $(s) and /3 = 0. 

- If on an open set, grad 9 A grad s 7^ then 

(2) =^> grad # A grad s = grad a A grad 
Consequently, grad#, grads, grad a and grad/3 are in the same plane and 
there exist tp, ip, x such that (3 = (p(9,s), a = ijj(9,s) and 7 = x($> s )- So, 

(*) («)■ 

(m) yields locally 9 = h(/3,s). 
There exists 7 = H(/3, s) such that H' s = 9. With a = —HL we obtain 
grad* j = H' s grad* s + H'^ grad* /3. So, (22) =>- (i). 

Lemma 2 Let 0~ o be an arbitrary non-constant differential scalar field in T> 
independent of a given scalar field 5 Q defined on T> by the relation <5 D (X) = 
%t(X)]. For any form field X G T> — > E*(X) G T£(£> ) i/iere eaisi -at 
least locally- three scalar fields defined on V Q and denoted by a a , j3 Q and uj 
such that 

S* = grad* a Q + f3 Q grad* 5 Q + u Q grad* a Q 
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Due to the fact that grad* S Q and grad* a D are two independent forms of V Q , 
the two surface-families S (X) = a\ and o~ Q (X) = b\ (where a\ and b\ are two 
arbitrary constants), generate a curve family denoted by (T) aibl . 
Let (S) be any transverse surface of the curve family (T) aibl , point M Q denotes 

the current intersection between (r) ai & 1 and (S). We denote by M Q M the arc 
of the curve (r) fll { )1 with origin at M Q . 
We define a differential function as 



XGD ^a (X)= f n E* dX 

JMoM 



such that d5 Q = 0, da Q = and da Q = S* <iX. Then, there exist -at least 
locally- two scalar field forms f3 Q and u in T> such that 

H* dX — da Q = (3 d5 Q + uj q da Q . 

Theorem 3 For any mapping o ofV t , independent of 5 and verifying — = ; 
for any form field H* in T> t with a zero convective derivation, there exist -at 

least locally- three mappings a, (3 and u ofV t , verifying = 0, ^ = and 

' dt dt 

—— = such that E* = grad* a + (3 grad* 5 + u grad* a. 
dt 

Let us consider a (t,x) = a Q [^ 1 (X)], a(t,x) = a Q [^ t -1 (X)] , /3(t,x) = /^[^(X)] 
and the field H* of forms in T> t corresponding to S*. 

Since a, a, (3, u and 5 have zero material derivatives, the previous forms and 
scalar functions have a zero convective derivatives. 

Consequences: 

From the representation of the form field C* and for 5 = T, we can get two 
scalar fields ip and ip in System (flOl) ; we obtain the potential representation 
(or Clebsch's representation) for the flows in System (Tl2|) . We note that the 
choice of C* is arbitrary and ip and ip are not unique. 

Proof of system (jTSj) : 

System (1181) is directly obtained when G is independent of T. In the converse, 
only a modification is needed: for a given form field 5* in V Q , there exist -at 
least locally- two scalar fields a and (3 Q such that 

2* = grad* a Q + (3 Q grad* T 

It is an analog representation as Clesch's [16]. In this case, we consider a 
transverse curve (r) of T-surfaces and a scalar field 

a (X)= f n ~* dX 

J MoM 

defined on T> as the circulation of H* on an arbitrary curve on the T-surface 
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connecting point X at point X G where the surface intersects (T). 
The circulation is independent of the curve connecting X G and X; when dT Q = 
then da = 2* dX. Consequently there exists a Lagrangian multiplier /3 (X) 
such that 

2* dX. = da Q + (3 dT Q 
and the Clebsch form is obtained. 
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